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Neutral silicon clusters Si TO (up to n = 7) and their cations SiJ (up to n = 10) have been 
produced by femtosecond laser ablation of bulk silicon in vacuum and investigated using time-of- 
flight mass spectrometry. Two populations of the Si^ clusters with different velocity and abundance 
distributions in the ablation plume have been clearly distinguished. Possible mechanisms of cluster 
formation (Coulomb explosion, gas-phase condensation, phase explosion) are discussed. 
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Femtosecond laser ablation is a rapidly developing 
technique offering new possibilities in various applica- 
tions. The fundamental mechanisms that lead to ma- 
terial removal are, however, still poorly understood. The 
interaction of fs laser pulses with silicon is an example 
of a complex interplay of thermal and ultrafast, non- 
thermal processes involved to ablation 0,0,13. Studies of 
composition and expansion dynamics of the laser-induced 
plume can provide a considerable insight into the ablation 
mechanisms. Several recent experiments on this subject 
0,0, IE El were mainly devoted to the atomic component 
of the plume though an observation of small silicon clus- 
ters was also mentioned ■ Little is known about silicon 
cluster formation under ablation with short laser pulses. 
With ns pulses, low-fluence desorption of Si dimers 
and small neutral Si„ clusters were observed. Desorp- 
tion of small clusters from both crystalline and nanos- 
tructured Si surfaces was induced by using high-energy 
(6.4 eV) photons 0. Multiple-charged cluster ions were 
formed in ps-laser stimulated field evaporation |lfj| . Re- 
cently we reported the first results on the observation 
of neutral and cationic clusters under femtosecond laser 
ablation of silicon ^lj. In this work, we analyze mech- 
anisms of clusters formation based on measurements of 
the abundance and velocity distributions. 



I. EXPERIMENT 

The experiments were performed with Si[100] surface 
under ultrahigh vacuum conditions (~ 10~ 10 mbar). The 
Si target was irradiated at a 45 ° incidence angle using 
a Tirsapphire laser (Mai-Tai coupled with a TSA ampli- 
fier, Spectra Physics, 80 fs pulse duration, 10 Hz rcp- 
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ctition rate, up to 30 mJ energy per pulse) operating 
at 800 nm. A part of the laser beam was selected by 
an aperture to provide a nearly uniform intensity dis- 
tribution over the irradiated spot. The target was ro- 
tated/translated during measurements to avoid consid- 
erable cratering. Some experiments were performed with 
the fixed target in order to investigate the effect of ac- 
cumulation of laser pulses at the same spot. The fluencc 
on the target was varied in the range 80-800 mJ/cm 2 . 

The abundance distributions of neutral and positively 
charged particles in the laser-ablation plume were stud- 
ied using a reflectron time-of-flight mass spectrometer. 
The neutral particles were analyzed using electron im- 
pact ionization (110 eV) and a plasma suppressor [T^ . 
At a distance of 11 cm from the target, the ablated or 
post-ionized ions were sampled parallel to the plume axis 
by a 500 V repeller pulse at a time delay td in respect 
to the laser pulse. Mass spectra were averaged over 300 
laser shots with the rotated target and over 5 shots with 
the fixed target. 



II. RESULTS AND DISCUSSION 

Monatomic Si + ions and Si atoms are the most abun- 
dant particles in the plume throughout the laser fluence 
range studied. The ion signal appears at a threshold flu- 
ence F th w 100 mJ/cm 2 [ll| for the etched Si surface 
irradiated previously with a fairly large number of laser 
pulses. Efficient emission of both neutral silicon clusters 
Si n (up to n = 7) and their cations Si+ (up to n = 10) 
has been observed for F > F th . The cluster abundance 
distributions are found to be essentially different for the 
fresh and etched ablated surfaces. Figure ^ shows a typ- 
ical evolution of cation mass spectra with increasing the 
number N of laser pulses applied to the same spot. The 
spectra were obtained at time delay td, corresponding 
to a maximum average yield of clusters. With the fresh 
surface, the relatively large Si+ clusters (n > 3) are ob- 
served in high abundance with the " magic" number at n 
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= 6 (Fig. [IJi,). In addition, doubly charged cluster ions 
are present in the plume among which Si| + is particu- 
larly abundant (note that Si 2 , -1 " was observed earlier as a 
magic cluster in the field evaporation experiment jlCjV 
The monatomic Si + peak is relatively weak under these 
conditions. 

When the number of laser shots increases, the mass 
spectrum is significantly modified (Fig. ^,0). The Si+ 
(n > 3) and Si§ + clusters are observed in lower abun- 
dance (the latter almost disappears at N > 100) while 
the concentration of the smaller species (n = f-3) in- 
creases progressively with N. After ~40 laser shots, Si + 
ion is totally dominant in the plume (this is not seen in 
Fig. n smce the Si + signal is maximized at shorter time 
delays). The cluster distribution exhibits an odd-even 
alternation with preferred formation of even-numbered 
Si+ clusters. Further increase in the number of laser 
pulses results in further decrease of the relative yield of 
the cationic clusters. After few hundreds of laser shots 
the signal saturates and the abundance distribution be- 
comes rather smooth with peak intensities monotonously 
decreasing with cluster size (Fig. QJ). 

Typical mass spectra for neutral Si„ clusters obtained 
with the fresh and etched ablated surfaces are shown in 
Fig. [21 The effect of the number of accumulated laser 
pulses is not so strong in this case though the same ten- 
dency of decreasing in the relative yield for Si„ clusters 
(n > 4) is obvious. The odd-even alternation in the clus- 
ter distribution is also observed at low N values (Fig. 
but it is less pronounced than for cationic clusters. 

The observed effect of the number of applied laser 
pulses on particle emission suggests an accumulation of 
the laser-induced damage on the Si surface. Indeed, with 
an optical microscope we observed clearly shaped spots at 
the surface for very low fluences down to approximately 
150 mJ/cm 2 (for multi-pulse conditions with ~I000 laser 
shots per spot). We can thus conclude that the dam- 
age threshold is approximately equal to the ion desorp- 
tion threshold for silicon under these conditions. The 
low-fluence damage (modification) of Si surfaces with ul- 
trashort laser pulses was observed earlier and attributed 
to amorphization of silicon [T3L ll4T | . The effect of laser 
pulse accumulation on the damage threshold (the incu- 
bation effect) was also revealed. The multi-pulse damage 
threshold determined in ^4| is well consistent with that 
evaluated here. Similar incubation effect in silicon was 
observed recently with ns laser pulses 0, ^| . The pre- 
cise nature of the incubation effect is not clear (it can 
be, e.g., due to chemical modification, mechanical stress, 
defect accumulation 0, 0] ) and needs further investi- 
gation. 

Variation of the time delay td allowed the analysis of 
particle velocity distributions and the characterization of 
the temporal evolution of the plume composition. The 
measurements were perfermed with the etched Si surface 
when no signal change with the number of laser pulses 
occurred. Figure shows the time-of-flight (TOF) dis- 
tributions of several cationic particles for different laser 
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FIG. 1: Mass spectra of cationic silicon clusters produced 
at laser fluence F = 350 mJ/cm 2 with different number N of 
laser pulses applied. Time delay td = 32 /is 



fluences. The distributions for Si + ions are maximized at 
^9 /is time delay that corresponds to ~12 km/s ion ve- 
locity. Increase in fluence results just in a broadening of 
the distributions with little or no shift towards higher ve- 
locities. At the same time, the yield of Si + ions increases 
strongly with fluence as ~ F n with n m 7.5 [Tlj . 

In contrast to Si + ions, the cluster TOF distributions 
are found to be strongly fluence-dependent. At low flu- 
ence, near the threshold value, the TOF distribution for 
S12 is rather narrow, single-peaked, and maximized at 
td = 18 /is (Fig. |3Ji), that is the most probable velocity 
is approximately two times lower than that of Si + . This 
means that these plume particles have nearly the same 



momentum. At 400 mJ/cm , the distribution for Si J is 
still single-peaked but a noticeable shift towards higher 
time delays is observed (the maximum occurs at ^25 
/xs, the corresponding velocity is ~ 4.4 km/s). The fast 
SiJ dimers are still present in this ablation regime but 
their distribution is masked by slower ions. With further 
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FIG. 2: Mass spectra of neutral silicon clusters produced at 
F = 350 mj/cm 2 and detected at ta = 32 fis with relatively 
fresh (a) and etched (b) surfaces. 



increase in laser fluence, the second, even slower, popu- 
lation of SiJ appears in the plume. At 550 mJ/cm 2 this 
slow population dominates (Fig. The first faster 

peak in the distribution is still present, well separated 
from the second peak, and is still maximized at ^25 fis. 

The evolution of TOF distributions for larger Si+ clus- 
ters is qualitatively similar to that for SiJ as illustrated 
in Fig. U] for the Sg" cation. At a threshold fluence of 
around 450 mJ/cm 2 , the distributions are transformed 
from single-peaked to double-peaked. The second (slow) 
cluster population becomes rapidly dominat with further 
increase in fluence. The most probable velocity of the sec- 
ond population decreases slightly with cluster size from 
2.5 km/s for SiJ to 2 km/s for Sig". The abundance 
distribution for this slow cluster population exhibits a 
pronounced odd-even alternation with dominance of the 
Si^~ peak and similar to that observed at shorter delay 
times with relatively fresh surface (Fig. QJi). It is thus 
remarkably different from the falling smooth distribution 
for faster clusters (Fig. [TJ;). 

Based on the obtained results, we can address the 
fundamental question on cluster formation mechanism 
under different irradiation regimes. At very low laser 
fluence, from the ion appearance threshold F t h up to 
^2F t h, Si + and Si J ions are likely produced by an im- 
pulsive Coulomb explosion (CE) from a charged surface. 
Two observations support the CE mechanism Il6j| . 
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FIG. 3: TOF distributions for SV 
three laser fluences. 
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namely, (a) the momentum scaling for the particles un- 
der low fluence regimes and (b) decrease of velocity for 
Si^" ions when fluence is increased beyond ^2F t h, that is 
a behavior opposite to thermal desorption. The slower 
Si+ clusters (which overshadow the CE ions in the (2- 
4) fluence range and form the fast cluster population 
at higher fluences) can be interpreted as "plasma ions", 
i.e., ions formed in the ionized vapor plume by collision- 
induced condensation. The velocities of these "plasma 
ions" decrease only slightly with cluster size and scale by 
a law intermediate between constant kinetic energy and 
constant velocity that implies the gas-phase condensation 
mechanism rather than direct ejection of the clusters [l2l |. 

Of particular interest is the origin of very slow clus- 
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ter ions observed at high fluences (F > A.hFtu) as the 
second population in the TOF distributions (Fig. 
The evident separation of this population from the first 
one indicate that another mechanism is likely respon- 
sible for its formation. We suggest that the slow ions 
are due to phase explosion (PE) from the Si surface 
melted via an ultrafast, non-thermal process which was 
recently found for semiconductors under fs-laser irradia- 
tion P, 0, 01 • The fast excitation of a dense electron-hole 
plasma causes destabilization of the lattice, and the semi- 
conductor melts in a time less than lps. The non-thermal 
melting occurs at fluences fairly higher than needed for 
purely thermal melting process. The melted surface, even 
if its temperature is far below the thermodynamic critical 
temperature, inevitably enters the region of metastable 
states with the subsequent phase exp losion as one of 
probable stabilization mechanisms |l7j . The overheated 
liquid can decay into the gas phase without formation 
of the critical vapor nucleus [l7j so the ejection of small 
clusters is quite possible. More work is needed to be 
certain about the origin of the cluster ions under these 
conditions but some observations support the PE mech- 
anism. The amount of material released under the PE 
is macroscopic and the ejected particles undergo many 
collisions within the plume during their expansion into 
vacuum. As a result, the most stable clusters survive 
predominantly in the plume thus resulting in abundance 
distributions with magic numbers. Numerous previous 
investigations indicate special stability of SiJ, Si^T , and 



particularly Si^j" ■ This is well correlated with the ob- 
served abundance distributions of the slow ions (similar 
to that shown in Fig. In addition, the observed 

weak size dependence of cluster velocities for the second 
population (Fig. |3J;) correlate with recent simulations 
of short-pulse laser ablation where nearly the same 
flow velocities were obtained for clusters of different size 
ejected via the PE mechanism. 



III. CONCLUSIONS 

We have observed efficient emission of both neutral and 
cationic silicon clusters under femtosecond laser ablation 
of silicon in vacuum. Cluster abundance and velocity dis- 
tributions are found to depend strongly on laser fluence 
and the number of laser pulses applied. The obtained 
results provide clear evidence that different cluster for- 
mation mechanisms (Coulomb explosion, gas-phase con- 
densation, phase explosion) are involved. Further work 
aimed at clarifying the contributions of these mechanisms 
for various ablation conditions is currently under way. 
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